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Abstract: We present a standard hydrodynamical description for non-canonical scalar 
field theories with kinetic gravity braiding. In particular, this picture applies to the simplest 
galileons and k- essence. The fluid variables not only have a clear physical meaning but 
also drastically simplify the analysis of the system. The fluid carries charges corresponding 
to shifts in field space. This shift-charge current contains a spatial part responsible for 
diffusion of the charges. Moreover, in the incompressible limit, the equation of motion 
becomes the standard diffusion equation. The fluid is indeed imperfect because the energy 
flows neither along the field gradient nor along the shift current. The fluid has zero vorticity 
and is not dissipative: there is no entropy production, the energy-momentum is exactly 
conserved, the temperature vanishes and there is no shear viscosity. Still, in an expansion 
around a perfect fluid one can identify terms which correct the pressure in the manner of 
bulk viscosity. We close by formulating the non-trivial conditions for the thermodynamic 
equilibrium of this imperfect fluid. 



Contents 

1 Introduction 1 

2 Action and Equations of Motion 5 

3 Fluid Picture 7 

3.1 Four- Velocity and Kinematical Decomposition 7 

3.2 Shift-Charge or Particle Current 8 

3.3 Chemical Potential and Force 9 

3.4 Energy-Momentum Tensor 11 

3.5 Energy-Momentum Conservation and the First Law of Thermodynamics 12 

3.6 Equation of Motion 13 

3.7 Diffusion and Non-Dissipative Bulk Viscosity 15 

4 Equilibrium and Vacuum Configurations 19 

5 How unnaturally complicated could it be? 20 

6 Discussion and Future Directions 22 
Acknowledgments 24 
A Choice of Frames 24 
B Action Without Second Time Derivatives 26 
References 27 



1 Introduction 

Scalar field theories with non-canonical derivative interactions, often referred to as essence 
[1, 2], have gained significant attention in the past decade as novel and phenomenologically 
different models for the inflationary stage of the early universe [3, 4] and of the recent 
universe, dominated by dark energy ("DE"), see [1, 2, 5] and dark matter ("DM"), see e.g. 
Refs [6-13]. The key new observable appearing in these theories in comparison to canonical 
scalar fields is the speed of sound which can be very different from the speed of light. 

On the other hand, modifications of gravity like [14, 15] in the so-called decoupling limit 
[16-18] ^ reduce to non-canonical scalar field theory structurally different from k- essence. 
In particular, in Minkowski spacetime this scalar field theory by construction possesses 

""^For more recent development, see [19, 20] 
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the galilean symmetry: — >■ d^(p + where = const. A (/ah'/ean-symmetric general- 
isation of this scalar field theory was introduced in [21, 22]^. Further, it was found that 
it is impossible to covariantise these models in such a way so as to maintain the galilean 
symmetry for dynamical spacetimes in a self-consistent manner [26, 27]. Nonetheless, these 
theories have been named galileon theories and a significant body of work has already been 
carried out analysing their impact not only on cosmological solutions [28-45], but also in 
other circumstances [46-55]. In [56], we introduced a class of models which extends both 
\l- essence and the action of the DGP [14] decoupling limit by 



which we have utilised as an Imperfect Dark Energy (see also Ref. [57], where this class 
of models was slightly later studied under the name G-inflation) This action does not 
possess the galilean symmetry and yet it does not contain a new degree of freedom, even 
in the presence of gravity. We have found a number of surprising features exhibited by 
this class of models. For example, contrary to k-essence and the ghost condensate [5], 
even in the limit where the scalar is an exact Goldstone boson, without direct couplings to 
matter, there exist attractors in expanding cosmologies in which the scalar field monitors 
and responds to the external energy density, only to eventually arrive at a final de-Sitter 
state; the null energy condition is generically violated in a stable manner and the system 
can evolve so as to cross the phantom divide. These features arise as a consequence of what 
we have named kinetic gravity braiding, the essential mixing of the derivatives of the scalar 
and of the metric, which cannot be undone by a field redefinition and which necessarily 
modifies gravity. The result of this kinetic braiding is that the energy-momentum tensor 
can no longer be brought to the perfect-fiuid form. A number of works have followed up 
the analysis, investigating structure formation and inflation in this class of models [58-62]. 

In this work we extend our heretofore investigation of the cosmological solutions to 
generic backgrounds and analyse the models with kinetic gravity braiding as fluids. This 
turns out to be a very fruitful framework in which to proceed and we demonstrate that 
it results in a very signiflcant simpliflcation of the equations and in giving a very physical 
picture to the a priori complex dynamics of the system. For the complete list of fluid 
variables and notation see table 1. 

We begin by recapping parts of our discussion on the equations of motion already 
presented in Ref. [56] in section 2. Then, in section 3, we show that a fluid description is 
still possible. Just as in the case of k- essence, one identifies the time-like gradient of the 
scalar field with the fluid's velocity in this way choosing a local rest frame. This implies 
that the scalar fleld plays the role of the internal clock for the fluid. Since the shift- 
symmetric Lagrangian is necessarily asymmetric with respect to — t- —(p there is a built-in 
arrow of time. Secondly one can realise that the equation of motion for the scalar takes 
the form of a divergence of a current, which is conserved in the case of a shift-symmetric 
Lagrangian. This allows us to identify the shift-charge density, n, in the local rest frame. 

^For earlier works unrelated to the decoupling limit, see e.g. [23-25] 
^For a further generalisation of galileons and our action (1.1) see Ref. [53] 
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What will prove key in the discussion here, is that this frame is neither the frame in which 
the charges are at rest, nor is it the frame in which there is no energy flux. However, it 
is in this particular frame that vorticity vanishes providing the hypersurface of constant 
intrinsic time on which the Cauchy data can be posed. 

We follow Schutz [63] and interpret the derivative of the scalar field with respect to 
proper time as the chemical potential, m. The surprising new feature is that both the 
charge density and the the energy density now explicitly contain the expansion of the fluid 
elements. In addition to its thermodynamical part P (m), the pressure, V, contains a term 
proportional to the first time derivative of the chemical potential. As we will show, these 
terms are a signature of the fluid's generically being out of thermodynamical equilibrium. 
In addition, this fluid is now imperfect. Indeed, the energy flows neither along the field 
gradient nor along the shift current.^ Therefore the energy-momentum tensor contains off- 
diagonal energy- flow terms. This energy flow is, in fact, diffusion occurring along gradients 
of the chemical potential. 



where V is the spatial gradient and k is the diffusivity. This is a dissipationless form of 
diffusion, with no entropy production. The fluid in this picture has zero temperature, is 
vorticity-free and does not have any shear viscosity, just one would expect of a system 
described by a single scalar field and an action. Thus it is natural to identify this fluid 
with an imperfect superfluid. The use of the fluid variables allows us to understand the 
dynamics described by the highly nonlinear Ampere-Monge-like equation of motion as the 
simple first law of thermodynamics 



where V is the comoving volume and M is the number of shift charges contained in V. 
In particular, the function K in the Lagrangian (1.1) provides the equation of state for 
the fluid, while the function G describes the dependence of the diffusivity on the chemical 
potential. Note that the usual Galileon choice for the function G = (dcj))'^ is not necessarily 
one which is physically motivated in this fluid picture, as we discuss in section 5. 

In section 3.7, we show that when the fluid dynamics are expanded around the perfect 
fluid, with the diffusivity k employed as the expansion parameter, one obtains a correction 
to the pressure that is proportional to the expansion of the fluid elements, and therefore 
behaves in the manner of bulk viscosity. However, when a proper gradient expansion is 
taken, we show that this bulk viscosity disappears, confirming that there is no dissipation. 
That being said, the fluid does respond to the expansion, and therefore we propose that 
it is viscid in this generalised, non-dissipative sense. Secondly, in the incompressible limit, 
the equation of motion for the scalar field can be re-expressed as the diffusion equation for 
charge density 



q = —K.'Vm 




dE 



VdV + mdhf , 



(1.3) 




(1.4) 



*Note the this definition of imperfection is different from the one used in Ref. [64, 65]. 
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with the diffusion coefficient T) = — k {mdn / dm)~ , where we have ignored the terms 
arising from the non-inertial nature of the rest frame. Yet again we reiterate: this diffusion 
process is dissipationless and occurs at zero temperature. 

Using the fluid language for \.- essence, we show that what is usually considered to 
be the relation between pressure and energy density (the "equation of state") is in fact 
nothing but the Euler relation for the thermodynamical representation of the fluid, relating 
the energy and pressure with the particle number and chemical potential. This relation 
guarantees the conservation of momentum. Surprisingly, the Euler relation is still valid for 
our fluid, but one must replace the total pressure with only its thermodynamical part. 



Table 1. Comparison of properties and notation of k- essence with the imperfect fluid of kinetic 
gravity braiding in the shift-symmetric case. 



Description 


Notation 




Definition 


Standard Kinetic Term 


X 






Fluid Velocity 






V„(I)/V2X 


Chemical Potential 


m 






Expansion 


e 




V„n^ 


Spatial Projector 








Time Derivative of x 


X 




nX 


Description 


Notation 


k- Essence 


Kinetic Gravity Braiding 


Lagrangian 


C 


K{X) 


K{X) + G (X) □(/. 


Diffusivity 


K 





2XGx 


Thermal Pressure 


P 


K 


K 


Charge Density 


n 


dP/dm 


dP/dm + kO 


Charge Current 


J, 


nu^ 


nUfj, — (k/tji) V yim 


Total Pressure 


V 


P 


P — Km 


Total Energy Density 


E 


nin — P 


mn — P 


Energy Flow 


% 





= -K _Lj; V^m 


Diffusion Coefficient 







—K (iTidn / dm)~^ 


Force Density 








Energy Conservation 




dE = -VdV 


dE = -VdV + mdN 
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In section 4, we turn to a discussion of the conditions required for thermodynamic 
equihbrium in our braided superfluid. A new result here (which also applies to k- essence, 
i.e. the usual perfect superfluid) is that, assuming our fluid interpretation, thermodynamic 
equilibrium implies the so-called Born-rigidity condition: the expansion and shear vanish; 
it turns out that the configuration must also be static. However, the key difference between 
these configurations and those considered previously (see e.g., [8]) is that the gradient of 
the scalar field remains time-like: therefore these type of equilibrium configurations only 
exist in theories with shift-symmetric Lagrangians. This also directly follows from our 
interpretation of (p as the intrinsic clock. 

A new feature arising from kinetic gravity braiding is that the equation of motion for the 
scalar field provides a new condition which must be fulfilled by equilibrium configurations. 
In particular, this condition constrains the equation of state which the scalar must adopt. 
In the extreme case of a constant diffusivity k, the total equation of state for the braided 
scalar together with external matter must he w = —1/3. In equilibrium, the scalar must 
arrange itself in such a way so as to screen the gravitational effects of itself and any external 
mass present in the system. 

We close by describing the motion of degenerate relativistic fermions by a k- essence 
scalar-field theory in 5. While the former is a simple physical system, the k- essence de- 
scribing it appears highly technically unnatural. 

In table 1, we are providing a quick dictionary defining our terms and contrasting the 
hydrodynamical properties of k- essence and the fiuid with kinetic gravity braiding. We 
assume the shift-symmetric case where the analogy with hydrodynamics is exact. 

2 Action and Equations of Motion 

In this section, we recap the discussion presented in our work [56, §2] where we introduced 
a class of scalar field theories minimally coupled to gravity, described by the action^ 



and K (</>, X) and G {<j), X) are arbitrary functions of the scalar field and its standard 
kinetic term X. Here and in most of the paper, we use reduced Planck units where 
Mpi = (SttGn)"^^^ = 1- For further discussion, it is convenient to use the standard 
notation for the corresponding Lagrangian 



which we consider as a function of three independent variables (j), X and B. One should note 
that, if the system is symmetric with respect to constant shifts in field space: (/>—>•(/) + c, 
then the Lagrangian is necessarily not symmetric with respect to cp ^ —(p. 

■^we use the metric signature convention (H — ). 





where we have denoted 
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The Lagrangian V equivalent to C can be obtained by integrating by parts the scalar- 
field contribution of the action Eq. (2.1): 



V = K 



G = K- 2XG4, - GxV^cpVxX , (2.3) 



where the subscripts (p and X denote partial differentiation with respect to these independ- 
ent variables. It is clear that it is the dependence of G on the field's gradient, X, that 
prevents this theory from being recast in a ]s.- essence form, whereas the derivative of G 
with respect to (j) contributes as an additional k- essence term. 

In Ref. [56], we showed that the equation of motion for the scalar field takes the 
following form 

J^' = , (2.4) 



where Va, = Kj, 



G(j, and the current 



= {Cx - 2G^) V^cj) - GxV^X , (2.5) 



corresponds to the Noether current for Lagrangians invariant under constant shifts in 
field space: (j) ^ (j) + c. Thus this current is conserved provided that = 0. The 
fully expanded equation of motion is somewhat unwieldy and was presented in Ref. [56], 
see Eq. (2.5) on page 6; it will not be necessary in the discussion here. We would like 
to stress that this second-order partial differential equation in the single scalar variable 
(f) is the only dynamical equation besides the Einstein equations. In particular, this is 
the equation which one has to solve in numerical analysis. The presence of third-order 
derivatives might have been expected owing to the d'Alembertian term in the action; 
however, as we have proven in [56], the equation of motion contains at most second-order 
derivatives. Therefore theories with kinetic gravity braiding do not contain any additional 
hidden (and ghosty) degrees of freedom which would appear by virtue of the Ostrogradsky 
procedure. Further, for any coordinate x*, the second derivative ViVicj) appears only 
linearly in the equation of motion. This allows one to associate this equation of motion with 
a particular generalisation of the Ampere-Monge equation for four-dimensional manifolds 
with Lorentzian signature. In section 3.6, we present a different derivation of the expanded 
equation of motion which makes use of the Raychaudhuri equation and explicitly shows 
how the higher-order derivatives are eliminated. 

The energy- momentum tensor ("EMT") for the scalar field is most easily derived in 
the standard way using the Lagrangian presented in Eq. (2.3) 

2 ^ S* 

T^, = = CxVf^(pV.^ - g,.vV - V,,GV,(t> - V,GV^<I) . (2.6) 

It is key to observe that the EMT also contains second derivatives of (j) appearing in 
"P, Cx and V^G. Therefore as we have already discussed in Ref. [56], both the Einstein's 
equations as well as the equation of motion for the scalar field (2.4) contain second deriv- 
atives of both the metric g^i^ and the scalar field (/>, so that the system is not diagonal in 
the second derivatives and no field redefinition exists which would unmix the two fields. It 
is this essential kinetic mixing of gravity and the scalar that we have named kinetic gravity 
braiding. 
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3 Fluid Picture 



When considering scalar-field theories as applied to cosmology, one usually restricts to 
the case with timelike field derivatives. It is well known that many aspects of k- essence 
in such a setup can be described in terms of relativistic hydrodynamics [3]. Below we 
introduce a fluid picture for the theory given by the action (2.1). As in the case of k- 
essence, this hydrodynamical/fluid language turns out to be helpful in developing intuition 
and in simplifying notation. To our best knowledge the analysis presented below appears 
in the literature for the first time. For earlier studies of thermodynamics of relativistic 
potential flows and k-essence see e.g. Refs [66-69]. 

3.1 Four- Velocity and Kinematical Decomposition 

First of all, by analogy with k-essence we can introduce an effective four- velocity 

^'"^ (3.1) 



^ V2X' 

which defines the local rest frame (LRF). We will then use the notation 

(•) = A( ) = n"V„( ), (3.2) 

for the derivative along (material derivative), thus making r the proper time of an 
observer comoving with the LRF. We can then introduce the transverse projector 

-L/ii/= Sfj.i' u^Ui, , (3.3) 

which will be the first fundamental form in the hypersurfaces : 4'{x) = const and which 
will allow us to decompose vectors (and gradients) into time-like and space-like parts as 
observed in the LRF, 

= u^u^Vx+ l-l Va = + . (3.4) 

Without loss of generality, we will assume that the time derivative of (j) is positive definite, 
always making future directed. The field (j) can thus be considered to be an internal 
clock. The hypersurfaces Ti^ then are Cauchy hypersurfaces for the standard-model fields 
and gravity and natural candidates for being the Cauchy hypersurfaces for the equation of 
motion of the field (f) (see the discussion regarding this particular choice of frame in appendix 
??). The equation of motion (2.4) is of the second order, therefore from the Cauchy- 
Kowalewski theorem, the initial data on are (f)Q {x) and (/jq (^c), which is represented not 
only by its value {4>q {x) = const on Tj^) but also by second derivatives with respect to the 
spatial coordinates a; on S^. 

In the shift-symmetric case, the Lagrangian is explicitly independent of the internal 
clock (j). Then, as we have already mentioned, the Lagrangian necessarily violates (j) ~^ ~4' 
symmetry furnishing our model with a built-in arrow of time or an asymmetry of the 
motions with respect to an increasing or decreasing internal clock (p. 
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The corresponding four-acceleration is defined as 



u^Vxu^ ■ (3.5) 



and can also be rewritten in form of a spatial gradient, 



o^=±^VAln0, (3.6) 



It will be convenient to make use of the notation for its magnitude, 



-a^ax . (3.7) 



Vf,Uy = u^a^ + -6 +a^i, , (3.8) 



Further, we can make the standard kinematical decomposition, 

1 

3* 

where 9 denotes the expansion, 

e Vxu'' = V^n'^ , (3.9) 



and aniy is the shear tensor. 



a^, = - Vxu,+ ±t VAn,.j - - ±^,u , (3.10) 

which is symmetric, traceless and purely spatial. In this decomposition, we have used the 
vanishing of the rotation tensor (or twist) VxUu— -L^ Vau^ in accordance with the 
Frobenius theorem. The four-acceleration a^, expansion 9 and the shear tensor cj^,y are 
constructed from spatial derivatives of the scalar field and can be calculated on any 
directly from the initial data without using the equation of motion. Finally, it is useful to 
introduce the extrinsic curvature of the hypersurface S,^, 

=-La-L^ VaU^ = 0"a/3 + ^ -La/3 , (3.11) 

from whence it follows that fC = 9. 

3.2 Shift-Charge or Particle Current 

Using the formulae 

X 



Vf,X = 2Xa^ + XUf, , and □(/> = ^= + V2X9 , (3.12) 

v 2>X 

we can write the Noether current Eq. (2.5), conserved only when shift symmetry is present, 
as 

= ((/) {Kx - 2G^) + k9) - , (3.13) 
where we have introduced the notation^ 

K = 2XGx ■ (3.14) 



®Note that in Ref. [56, Eq. (3.6)], we used a definition for k that is smaller than this one by a factor of 
2, Khcro = 2Kthcrc. This Served to simplify the equations even more in the particular case of a cosmological 
background, but the proper normalisation for this term which gives its physical meaning is the one presented 
in this work. 
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From the standard discussion of imperfect fluids (see e.g. Ref. [70]), we have to identify the 
space-like part of the current (3.13), —na^, with diffusion, while k can be interpreted as the 
diffusivity, which we will take to be the coefficient relating the energy flow to the gradient 
of the chemical potential'' (see Eq. (3.35)). Despite the presence of second derivatives in 
(2.5), the current does not contain the second time derivative of the field (p and can be 
calculated directly from initial data. The density of the shift-symmetry charge*^ 



n 



uf'J^ = HKx- 2G^) + , (3.15) 



contains the expansion 9, which is highly non-standard. It is convenient to associate this 
Noether current with the particle-number current. In such a way, the introduced total 
particle/charge number will be conserved provided there be shift symmetry. However, 
the number of particles M = nV in an infinitesimal volume V moving with velocity n'^ 
(comoving volume) is not conserved. Indeed using the decomposition (3.13) one can write 
equation of motion (2.4) in the form 

n + 9n- {Ka^') = , (3.16) 

and using the geometric meaning of the expansion 



(3.17) 



we can obtain 

^f = ViV^ + V^ (Ka'^)) = VV^ + Miif , (3.18) 

where A/dif denotes the number of particles transported to the volume by diffusion per unit 
of proper time. Thus generically Af is not conserved even in the shift-symmetric case. The 
total particle current is subluminal provided 



n 



^-K^a'^>0, (3.19) 



If Eq. (3.19) holds, then one can chose an alternative fiuid description where one takes the 
Eckart frame which moves with the particles as the LRF. We will comment more on this 
later. 



3.3 Chemical Potential and Force 

Using the general expression for the Energy-Momentum Tensor ("EMT") (2.6) and the 
decomposition (3.8), we define the energy density of the fluid in the standard way 

£ = T^^u'^u" = 2X {Kx - G^) + Ke^-K. (3.20) 

^Note that usually the terms diffusivity and diffusion coefficient are used interchangeably, since the same 
coefficient usually appears in both the Pick's laws. We have kept the definition of diffusivity as the one here 
and use it to mean a coefficient which determines an energy flux rather than a charge flux: in our case, this 
is a more natural definition. In section 3.7 we show that, in an expansion around an incompressible fluid, 
the equation of motion reduces to the diffusion equation with the diffusion coefficient related to, but not 
the same as, the diffusivity, see Eq. (3.59). 

*Note that in Ref. [56, Eq. (3.2)], we denoted this number density by J: Whcro = J'thore- 
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Contrary to the usual descriptions of relativistic imperfect fluids, the energy density con- 
tains the expansion, 6. 

For the shift-symmetric case we recall the definition of the effective mass per charge 
(which is also the chemical potential or injection energy) see e.g. [71, p. 562] and obtain 

™.(|£) =(|£) ,3.21) 

\Onj T-,y=const \Onj 0=const,6»=O 

this expression is identical to that of Ref. [63]. Further, we will identify 

m = = V2X>O, (3.22) 

with an effective rest mass per particle building the fluid even in non-shift-symmetric cases. 
As we will see later (see Eqs (3.37), (3.41)) the rest-mass introduced in this way plays the 
role of the chemical potential.^ 

The effective momentum of the particle can then be defined as 

= mu^, (3.23) 

and the corresponding relativistic vorticity, V ^^Pu — ^uP^, is zero.^'' Further, the fact that 
Pfj-P^ = and = V^cp, means that —cp plays the role of Hamilton's principal function 
(action as a function of final coordinates) for each shift-charge. Indeed, for the action of a 
particle we have 

mdT = -<j), (3.24) 



where we have used Eq. (3.22) and where we have omitted the constant of integration. 

Exploiting the particle momentum, one can introduce a relativistic 3-force acting 
on each particle. 



d 

^'f^ ^-l;: px =-l; 

Then, using the decomposition 



PX =±f, ^ {mux) . (3.25) 



V^m = ma^ + ?nn^ , (3.26) 

which follows from the definition of the effective mass (3.22), and decomposition (3.8), the 
force can be written as the gradient of the chemical potential: 

Ff, = ma,, Vxm . (3.27) 

Here it is worth remembering that, in a system in thermodynamical equilibrium which 
experiences an external potential, the chemical potential is given by this external potential 



^In the literature, m is also called specific enthalpy or Gibbs energy per particle. 

^°Here it is worth recalling that in a stationary gravitational field with Killing vector dt the conserved 
energy of a particle is governed by pt- Using our identification we obtain that pt — dtij) has to be conserved. 
As it was shown for k-essence in [72], the energy pt given by (3.23) is indeed conserved for stationary 
configurations. Taking into account that (j> plays the role of an internal clock one can expect that stationarity 
implies shift symmetry in field space. Thus we conjecture that the statement from [72] is also applicable to 
the kinetic gravity braiding theories. 
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plus some constant, so that the force acting on the particle is exactly the gradient of 
the external potential, as expected (do not forget the signature). We discuss equilibrium 
configurations in detail in section 4. 

Further we can also express the current Eq. (3.13) making use of the chemical potential 

J^, = nu^-- ±^ Vxm , (3.28) 
c.f. with the formula (59.5) for the diffusive current from Ref. [66, p. 231]. 
3.4 Energy-Momentum Tensor 

In this fluid language the EMT Eq. (2.6) can be re-expressed as 

T^u = £u^Uy- V + Ufj,qu + UuQf, , (3.29) 

where the energy density is 

£ = T^.uU^'u'' = 2X {Kx - G^) + Omn - K . (3.30) 

As we have already remarked, contrary to the usual descriptions of relativistic imperfect 
fluids, the energy density contains the expansion, 9. Still, this expression can be calculated 
purely from the initial data. 

Further, we can define the effective total isotropic pressure in the usual fashion and 
find that it is given by the Lagrangian Eq. (2.3): 

p = _ Ita*^ ^^y=K- 2XG^ -Km. (3.31) 

The presence of the last term containing m can be associated with an additional force 
transporting particles along the gradients of the chemical potential. Because of this term, 
the total pressure can be calculated only after making use of the equations of motion which 
is not the case in standard fiuids (see section 3.6 for a further discussion). However, in 
configurations in thermodynamical equilibrium this term vanishes. It will prove helpful to 
also introduce the part of pressure 

P = K- 2XG^ , (3.32) 

which can be directly calculated from the initial data and which does not vanish in equi- 
librium configurations, usually called the thermodynamic pressure 
Finally the energy fiow associated with diffusion is 

=-L^A T^u'' = -mna^ = m ±^i, J'^ . (3.33) 

Thus the energy fiow arises purely from diffusive part of the particle current _L^,y J"^ 
and the heat fiow, which is defined as""^^ 

Heat Flow =_L^a T^u" - m ±f,y J" = , (3.34) 

^^This imperfection from the construction called Inhomogeneous Equation of State where pressure is 
postulated to be a function of both the energy density and the Hubble parameter, see e.g. [73-75]. 

^■^See e.g. Ref. [76, §1, Eq. (39)]. Note that many works use the Eckart frame as the LRF from the very 
beginning so that the (/^ used there can only be the heat flow. 
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vanishes, as it should for a system with identically zero temperature. Then we can see that 
the diffusive energy flux is proportional to the spatial gradient of the chemical potential: 

q^ = -Kiy^Vum, (3.35) 

c.f. with the formula (59.5) for the diffusive energy flux from [66, p. 231]. We take this 
relation as our definition of k as the diffusivity. It is instructive to compare these relations 
with formulae (285) and (287) from the imperfect fluid from review [70, p. 58]. 

In simple equilibrium thermodynamics, the chemical potential has to be equal through- 
out the body. If there are gradients of the chemical potential, then the system develops 
forces and transports particles to compensate for these gradients. As we discuss in section 
4, in our case it is possible to construct equilibrium conflgurations which do nonetheless 
exhibit diffusive flow in a time-invariant fashion. 

Having deflned the thermodynamic quantities, we can then re-express the particle 
density Eq. (3.15) 

n ((/., m, e) = Pm + K9 + niK^ = V,n + (ku^) , (3.36) 

where we have explicitly stressed that we are treating the expansion as an independent 
variable. This form allows us to immediately see that n can be found using solely initial 
data. 

Using equations (3.15), (3.22), (3.30) and (3.32) we can obtain the Euler relation: 

S = mn-P . (3.37) 

Observe that the diffusion contribution nra cancels out and only the usual equilibrium 
part of the pressure, which can be calculated from the initial data, contributes. We should 
stress here that all the above results hold even when shift-symmetry is broken and particle 
number is not conserved. 

We would like to mention that the result Eq. (3.37) in the k- essence limit {Gx = 0) 
reduces to the formula sometimes referred to as the ]i- essence equation of state, £ = 
2XPx — P [1] . To our best knowledge, the relationship between this equation of state and 
the Euler relation appears to have been missed in the literature. 

3.5 Energy-Momentum Conservation and the First Law of Thermodynamics 

It is instructive to analyse the equations resulting from the conservation of the EMT 
V^Tlf = 0. For energy conservation we obtain: 

u^V^T^"" = £ + e{£ + V) - Vx (m/ca^) + rriKaxa^ = , (3.38) 

where we have used the definition of the expansion and the property of four acceleration 
u^a^ = 0. Further we can express the divergence of the diffusion current through the 
equation of motion (3.16) 

(ko^) = m (n + 9n - V^) , (3.39) 
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so that the conservation of energy eventually reduces to 

£ + 9 {£ + V) - in{h + 9n) + mV^ = . (3.40) 

Let us now consider the evolution of the energy E = £V in the comoving volume V. 
Recalling the definition of the expansion (3.17) and the "particle conservation" equation 
(3.18) we obtain the first law of thermodynamics for the matter in the comoving volume 
in the shift-symmetric case: 

dE = -VdV + m(W , (3.41) 

where TV is number of particles in the volume V, while d denotes the differential along 
the velocity u^. From expression (3.41) it is obvious that m plays the role of chemical 
potential. In the general case with "P^ ^ 0, one should modify Eq. (3.41) by substituting 
dM —7- SMdii — the infinitesimal number of particles transported to the volume by diffusion, 
but not those produced by the source V^f,. 

There is only one single scalar degree of freedom (p and one independent scalar equation 
of motion (2.4). Hence _L^jy \/\T^'^ = 0, momentum conservation, cannot contain any 
additional information. Indeed, one can calculate that 

^,.u ^xT^" = {£ + P-mn)a^, (3.42) 

therefore the Euler relation (3.37) guarantees momentum conservation. 

For completeness, it is worth mentioning how the "Archimedes law" changes for this 
fluid. The force density is given by 

= nF^ , (3.43) 
or applying the Euler relation (3.37) it can be written as 

= n ±^ Vatti = nma^ = {£ + P) . (3.44) 

On the other hand, using momentum conservation -L^^y V xT'^'^ = 0, after some algebra we 
obtain 

ff,=±^^VxP+{^^m + Ke) ±^^Vxm. (3.45) 

Hence the force density originates from a superposition of spatial gradients of the equilib- 
rium part of pressure and of the chemical potential. Using Eq. (3.36), we can rewrite the 
contents of the parentheses as K^m + k9 = n — Pm, which disappears in thermodynamical 
equilibrium (it is one of the Maxwell relations). On the other hand, in k- essence, the last 
term is not present by the simple virtue of k = 0. 

3.6 Equation of Motion 

Let us discuss the equation of motion (3.16) in more detail. In particular, we will be 
interested in solving it with respect to highest time derivative: m = (j). Naively it seems 
that Eq. (3.16) should have third derivatives. However, as we know from Ref. [56] these 
terms with third derivatives cancel. Indeed, using Eq. (3.26), (3.31) and the Euler integral 
(3.37) we obtain for Eq. (3.16) 

(rem + «^0) m + k{9 - V^a^j + 9n - niKma^ax + £^ = ^ , (3.46) 
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while from the Raychaudhuri equation we have 



e = -^9^ + V^a^ - a^^af'" - R^.u^'u'' , (3.47) 
so that third derivatives disappear giving a second-order equation 

{n„^ + K^) m - K (^o-„/30-"'^ + + R^yU^u"^ +en- K„,ma^ax + = . (3.48) 
We then use the Gauss equation 

cTa/3fT°'^ - = -^^^R - 2G^yU^'v!' (3.49) 

to re-express the shear in terms of the intrinsic curvature of the spatial hypersurface and 
obtain for the equation of motion (3.48) 

{rira + K^)m + dn + K {^'^^R + 2Gf,uU^u'' - Rf^uU^'u" - 6*^) - Kmma^ax + £4, = ^. (3.50) 

This equation of motion is linear in the highest time derivative (p = rh. Using Einstein 
equations and the expression for the EMT (3.29), we can eliminate the curvature term 

3 
2 



2G^^M^u'' - Rf.uu'^u^ = -{£-P + Km) + W^^t , (3.51) 



where we have denoted any possible contribution to the above contraction of the curvature 
tensors external to the scalar field as 

Wet = (tZ + ^T.^tg^'^) u^u, . (3.52) 



Thus after such an elimination of the second derivatives of the metric, the equation of 
motion takes the form 

nm + K,<t> + ^K^^ m+On+K (^^^ R + ^{£ - P) - 6*^^ - Kmma^afj,+£^ = -nW^t ■ (3.53) 

The diffusivity k also serves as an effective coupling constant in front of the external source 
term VFcxt- Since the total pressure V contains m (see Eq. (3.31)), the above equation 
implies that the pressure of the fluid depends on the acceleration and the extrinsic curvature 
tensor (and therefore the expansion 6). Moreover, we should expect that the coefficient in 
front of (j) is responsible for the presence or absence of ghosts / negative kinetic energies. 
Namely, there are no ghosts provided 

3 

D = n^ + K^ + -K^>^. (3.54) 

We confirm this claim explicitly in our discussion of the effective acoustic metric for per- 
turbations in the follow-up work Ref. [77]. As we have already mentioned, this condition 
is not correlated with the Null Energy Condition. Using Eq. (3.37) and Eq. (3.36) this 
function D can also be written in the form 

D = ^^li^ + Ik\ (3.55) 
m I 
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which simply reduces to the well-known 'k-essence limit, D = Ex [78]. Provided D 7^ 0, 
the solution of the Cauchy problem exists locally under the usual limitations. On the 
other hand, from the equation of motion (3.53) it follows that the zeroes of D generically 
correspond to pole-like singularities in m. Given the dependence of the total pressure on 
m, Eq. (3.31), the configurations with D = are pressure singularities and they build a 
barrier in phase space which is generically impenetrable. Thus a fluid evolving from some 
initial data with healthy perturbations (D > 0) will never become ghosty in a smooth 
and controllable way. Note that, similarly to the cases studied in Ref. [78], for some 
Lagrangians and some external matter Wext there may well exist trajectories penetrating 
this barrier. However, these trajectories would have measure zero in the phase space of the 
system since they correspond to an exact cancellation of the above singularity. 

Finally, we would like to stress that the fluid variables are helpful to develop intuition 
but do not appear to be very practical in actual numerical calculations. The equation of 
motion cannot be reduced to system of ordinary differential equations along a world line. In 
particular, the system (3.53) and (3.47) is not closed and requires the somewhat unwieldy 
Raychaudhuri equation for the shear tensor afj_u which then has to be supplemented with an 
even more complicated equation for d^. On top of that all these equations will be subject 
to the constraint that the vorticity of vanish. Working with (p directly and solving the 
expanded (2.4) has the advantage that there are no further constraints. Solving equations 
of motion for the velocity potential (phase for superfluids) (j) is of course the standard 
practice in vorticity- free hydrodynamics. 

3.7 Diffusion and Non-Dissipative Bulk Viscosity 

Let us now elaborate on the diffusive interpretation of this fluid. First of all, one quickly 
notices that Eq. (3.28) is a "1st Fick's law" J oc Vm, relating the diffusive flux of particles 
in or out of a comoving volume to the spatial gradient of the chemical potential, which 
again is consistent with our interpretation of k as the diffusivity. Of course, the "2nd Fick's 
law", namely a diffusion-type equation for the charge density n, can be identified with the 
equation of motion. Turning the attention to the shift-symmetric case for simplicity, we 
can rewrite the equation of motion (3.39) as 

n + en = rimm + 9n + Ke = V^ ±'"' Vurri^ , (3.56) 

which is already a time-dependent diffusion-type equation for m, with some extra terms. 
As we have already mentioned, the k6 term plays a crucial role to keep this equation of 
second order, through the Raychaudhury equation. One could worry that at the same time 
this may remove the diffusive character of Eq. (3.56), but this is actually not the case. For 
instance, there are configurations in which one can have a nonzero V^a^ but still ^ = (in 
Section 4 we discuss examples of such configurations). 

It is perhaps more illuminating to write Eq. (3.56) directly in terms of n. For this 
purpose, it suffices to realise that since n = n{m, 6) then at least locally one can express 
m = m(n,9). Performing a gradient expansion (i.e., in the expansion in n) one then has 
Vfj,m = :^Vi_in + 0{k) with n and taken at ^ = 0, that is for the perfect fluid limit. 
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With this, the equation of motion (or charge conservation) (3.16) reduces to 



h + en = Vx[ ^'"^ V^?i + O(k'). (3.57) 

Bar the terms involving the expansion, this equation is the standard diffusion equation 

n = -V^ (S) V'^n) + Tiaf^V^n , (3.58) 

with the diffusion coefficient a function of the diffusivity^^ 

^ = -^ = --c!^-T^cl (3.59) 
rimm n 

with Ts) the diffusion timescale appearing here since the diffusion coefficient has dimensions 
of L^/T. In the hmit of an incompressible fluid, Eq. (3.58) is in fact the exact equation 
of motion to all orders in k. The last term in this diffusion equation is proportional to 
acceleration and manifests that the equation is written in a non-inertial accelerated frame. 

Note that in this discussion we have not used the Einstein equations. The whole 
procedure is based on the presence of a higher derivative in the equation of motion, and 
the corresponding V^a^ in the equation of motion i.e. charge conservation (3.16). 

One may wonder whether the system gives enough freedom to ever realise this diffusive 
behaviour. It can be seen that this can occur when the initial data are such that the 
expansion vanishes and 6 is tuned to be small. This means that the timescale for the 
expansion of the volume is much larger than the timescale for the diffusion, 

rj, < . (3.60) 

This condition also happens to imply that the contribution from the expansion to the 
particle number is negligible. Then we must require that, in the charge-conservation equa- 
tion (3.16), the expansion term is negligible compared to divergence term, or roughly 

On < KV/,a^ ^ V^a^ > r^^ . (3.61) 

Finally, one needs to make sure that the configuration is maintained for sufficient time, i.e. 
that evolves slowly enough. Through the Raychaudhuri equation (3.47) we require that 

n 



(3.62) 



where all the above quantities depend purely on initial data. The final condition is pre- 
cisely the tuning that we are required to perform on these initial conditions to ensure that 
diffusion occurs for some finite time much longer than the characteristic diffusion timescale 



^^It is useful to compare this expression for X) with [66, Eq. (59.9), p. 232]. Note that our ti/m — —a 
and pc = n. 
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This is a rather fine-tuned situation, which is already clear from the fact that this can- 
cellation is not possible in 1+1 dimensions. However, initial conditions obeying (3.62) 
definitely exist. For example, in 1+3 dimensional Minkowski spacetime, these initial 
conditions can be easily realised by the following ansatz in Cartesian coordinates (t,x): 
c)f0o(x) = 1 and (/>o(x) = e/(x), where /(x) is a harmonic function, A/(x) = 0, and e is a 
small parameter. In that case the geometric invariants are 6 = 0{e^), V^a^ = 0{e^) and 
(^'^'^cTfii/ = C'(e^)i whereas from the Raychaudhuri equation we get 9 = 0{e'^). 
We conclude that the interpretation of this system as a diffusive one seems appropri- 
ate and robust. We emphasise that this is a non-dissipative form of diffusion present 
at zero temperature. Further one can see that Sa'^V^n = S (^nj / (iniirn) where 

^Vn^ = — J-^^'^ ^f^n'Vun. So that Eq. (3.58) is in fact a nonlinear diffusion equation. 

Let us now turn to the identification of the imperfect part of the EMT, which can 
be done by comparing with the conventional gradient expansion performed in relativistic 
thermo/hydrodynamics. As we have already mentioned in our discussion of the EMT in 
section 3.4, it is necessary to solve the equation of motion for the scalar field in order to 
express the EMT purely through variables on the Cauchy hypersurface. In doing so, we 
rewrite the EMT in the effective form of a gradient expansion with k, playing the role of 
the expansion parameter. Substituting the equation of motion (3.53) and the expression 
for the particle current in fluid variables (3.36) into the expression for the pressure (3.31) 
gives an EMT of the form 

T^,y = £u^u,y ~ ~ -^Mi' + '^''1^J, + C('t^) , (3.63) 

where we now have a pressure term corrected for the first-oder terms in this expansion, 

P^P+^£^, (3.64) 

and there appears a bulk viscosity-like term with a would-be viscosity coefficient 

C = -D-^KPm + 0{k^) = -K (—\ + 0{k^) , (3.65) 

\nmJ 

and we have neglected the terms second-order in the diffusivity k. We can compare this 
expression with the sound speed that we have obtained in [56, Eq. 3.20] for the cosmological 
background. Again dropping all terms of order k^, we can conclude that when the 
correction from kinetic braiding is small, the effective bulk viscosity, 

C = -Kmcl +0{k^). (3.66) 

However, it is important to bear in mind that this type of bulk viscosity differs from the 
usual notion defined in terms of the conventional gradient expansion. Indeed, in this expan- 
sion an important assumption is that the energy density does not receive any corrections 
involving gradients of potentials or velocities. In our fiuid, however, we do have a con- 
tribution to S- proportional to 9. Hence, the only way to appropriately compare to the 

^"^As we show in our forthcoming work Ref. [77], the sound speed on general backgrounds is more com- 
phcated but gives the same results at the level of precision required in this discussion. 
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language of the gradient expansion is to perform a redefinition of what is meant by the 
hydrodynamic potentials, order by order in k. Thus, one has to introduce a new chemical 
potential 

7h = m{m,9), (3.67) 

in such a way as to keep the form that is the energy density unchanged between the zeroth 
and first orders in the expansion, 

£o{m)=£{m,e) = £o{m) + Kme, (3.68) 

where Sq = mPm — P is the k- essence energy density. All quantities are supposed to be 
expressed now in terms of m. At leading order, we therefore have 

m K, m + — — = m + — . (3.69) 

the charge density n becomes a function of m only, up to first order in k 

n{m) = no(m) + kO = no{m) + 0{k^) . (3.70) 
This allows us to identify this modified chemical potential with the definition 

m = ^ , (3.71) 
on 

which is valid to the same order up to same order without the need to specify that the 
derivative is taken at constant volume. The pressure using the new potential becomes 

P(m) = P{m) + Kmcld . (3.72) 

When combined with the bulk viscosity term we obtained in Eq. (3.66), we can see that 
the contributions proportional to the expansion cancel out, and therefore the correct 
identification for the bulk viscosity of this fluid is 

Cgrad cxp — 0. (3.73) 

This is in tune with the fact that no actual dissipation takes place, as indeed it must not 
in the case of a theory arising from an action principle. 

Having said that, we must emphasise that the EMT does depend on the expansion 
9: even the energy density does (in the original variables). This corresponds to a kind 
of "viscous" behaviour: a (perhaps new) type of nontrivial response to expansion that is 
not dissipative. Hence, it still seems appropriate to call this a viscid fluid, even if the 
conventionally defined linear bulk and shear viscosity coefficients vanish. 

Taking as a measure of this type of viscosity the coefficient defined above in Eq. (3.66), 
we can see that there is also a relation between viscosity and diffusion in this fluid: 

C = 2)nm. (3.74) 



^^The same redefinitions of the chemical potential (3.67), particle number (3.70), energy density (3.1 
and pressure (3.72) were recently rederived in the work [79], c.f. equations (A9)-(A12). 
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4 Equilibrium and Vacuum Configurations 



In this section, we will discuss the conditions on the fluid that are implied by the require- 
ment that it is in equilibrium. As we will show here, the coupling of the braided fluid to 
gravity imposes an additional condition on the braided fluid compared to k- essence. 

Let us assume a non-vanishing number of particles in the fluid, n 7^ 0. In equilibrium 
in a stationary gravitational fleld with a timelike Killing vector the quantity which is 
constant throughout a system (in our case — on the spacelike hypersurfaces S^) is not the 
locally measured chemical potential m but 

m^Jgi^u^t^e = const. (4.1) 

This is sometimes called Tolman-Klein condition, see Ref. [80] and Ref. [81, §27, p. 77]. 

In the discussion here we will concentrate on the shift-symmetric case. Firstly, since (j) is 
an intrinsic clock and in equilibrium no physical observables should depend on this intrinsic 
clock explicitly. Secondly, comoving observables like £, m, etc. should be independent of 
comoving time, e.g. £ = u^V\£ = 0. Since we are building our description on the 
assumption that m = 7^ 0, it better be the case that (j) is not an observable, at least in 
equilibrium. 

The requirement that observables be independent of comoving time implies that the 
Killing vector .^^ oc u^. Since is hypersurface orthogonal by definition, this implies in 
fact that the equilibrium configuration is not just stationary, but also static. To satisfy Eq. 
(4.1) in an arbitrary curved spacetime automatically — without new conditions — we have 
to assume that the coefficient of proportionality is so that 

in equilibrium £,u = ~ is a Killing vector. (4.2) 

m 

Therefore 

V^^M + V,e. = ^ - = , (4.3) 

a projection of which gives 

m = and /C^i. = . (4.4) 

In particular, the expansion must vanish, 6 = /C^ = 0. Thus we have shown that, in 
equilibrium, the motion of the fiuid should be rigid in the Born sense: /C^,y = 0. This is 
the standard requirement, see e.g. Ref. [82]. 

We can also consider situations where more than one type of matter is present. Since 
^'^ is the only time- like Killing vector, the velocity of the external fiuid must lie along 
and since the velocities are normalised — tig^t = u^- The Tolman-Klein condition Eq. (4.1) 
then implies chemical potentials of all the constituents must be equal, ruext = fn- 

As we have already mentioned, the chemical potential m for the system in an external 
field can be space dependent even in equilibrium. Thus the acceleration does not necessarily 
vanish. However, this is a result of the fact that a non-zero acceleration implies that u^^ 
is not tangent to a geodesic: the energy flux in equilibrium appears since the LRF is non- 
inertial. One should stress, however, that the time dependence of the acceleration does 
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vanish. The appropriate measure of the time dependence of a vector in the Lie derivative, 
and in this, case 

£^a^, = a^V^ ( — ) + — Vac^ = . (4.5) 

In general, to find the equihbrium configuration as a function on space, one would need 
to apply the equilibrium conditions Eq. (4.4) to the Raychaudhuri equations and Einstein 
equations and solve this system. This is potentially a rich problem and lies outside of 
the scope of this work. Nonetheless, let us elucidate one significant difference between 
k-essence and the braided fluid. Under the equilibrium conditions, the equation of motion 
for k- essence is trivial and does not supply any new information. However, in the case of 
the braided fluid, the equation of motion (3.46) reduces to a new constraint 

kR^vU^v!^ = —Kmrna^ , (4-6) 

which, using the Einstein equations, becomes 

S + 3P + fext + 3Pext = -2(ln K)'a2 , (4.7) 

where the prime denotes differentiation with respect to Inm and the quantities with the 
subscript external are any potential contributions to the EMT coming from matter external 
to the braided scalar. Let us just say that this constraint imposes a condition on the 
equation of state of the braided scalar in equilibrium. In the extreme case of a constant 
diffusivity k, this constraint implies that the equation of state is such that the scalar 
arranges itself in a conflguration which screens the gravitational effect of the mass of 
the total EMT. In particular, when there is no external matter, the equation of state in 
equilibrium is w = —1/3. This is another facet of the monitoring behaviour present in the 
cosmological solutions which we discussed in detail in Ref. [56]. 

Let us now briefly mention vacuum configurations where the charge density vanishes, 
n = 0. In Ref. [56] we showed that for the shift-symmetric theories with kinetic gravity 
braiding there exist such non-trivial vacuum solutions and that they are dynamical at- 
tractors in an expanding universe. As one can see from Eq. (3.15), the vanishing of the 
shift-charge implies that the motion is restricted to obey 

0{m) = -—. (4.8) 

K 

The structure of these configurations is potentially even richer than that of equilibrium 
configurations and this issue deserves a further detailed analysis going far beyond of the 
scope of this paper. 

5 How unnaturally complicated could it be? 

One really interesting example of how unexpectedly complicated and unnatural the function 
K (X) can be is provided by the noninteracting relativistic degenerate fermions with mass 
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nie and chemical potential m, see e.g. [81, §61, p. 180]. In equilibrium, we have 
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where ^ is related to Fermi momentum pp as 



^ = 4arcsinh^^ . 
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(5.2) 



First of all, one can check that the Euler relation, Eq. (3.37), indeed holds. Further, we 
can re-express the pressure through the chemical potential m so that we obtain 
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If we restrict our attention to the simplest motion — vorticity-free — we have to follow 
Schutz [63] and the identification used in this work, Eq. (3.22), and substitute m = 

V2X = 



{d(f)) into this formula for the pressure to obtain the Lagrangian. As the result 

■esser 
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of this procedure, we get a rather unusual scalar field theory of the k- essence type: 
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Note that this Lagrangian describes a simple physical system and the incorporation of 
any interaction would only serve to make the structure more complicated. Moreover, 
this Lagrangian does not have a well defined vacuum limit at {dcf))^ — )• 0. This is the 
manifestation of the simple physical fact that the fermion chemical potential can never be 

n2 . „ 9 



smaller than the mass of the fermion me. Even in the limit {d4)Y 
not analytic: 



P {94') ^ TT^'TT-c 
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TTig the asymptotic is 



(5.5) 



Observe that the fundamental degrees of freedom in this system are fermions which have 
nothing to do with the effective bosonic field (p. Moreover, the Lagrangian does not appear 
to be technically natural and the only free parameter m^ does not correspond to the actual 
strong-coupling scale. One could expect that the latter is governed by the scale of particle 
separation, n^/^, which can be exponentially larger than me- 
lt is important to stress here that the Lagrangian Eq. (5.5) does not model all the 
degrees of freedom which exist in an ideal Fermi gas. In particular, it does not contain 
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the zero sound, which is a non-equihbrium propagating mode. In ideal Fermi gases this 
mode dominates at low temperature, while the adiabatic hydrodynamic mode modeled 
by Eq. (5.5) is actually suppressed. However, the addition of even a small attractive 
interaction allows the adiabatic mode to propagate [83, p. 163]. Adding such a small 
attractive interaction should have only a small effect on the equation of state of the fluid, 
since it only affects those modes lying close to the Fermi surface. On the other hand, we do 
not know whether there exist real examples of fluids where the hydrodynamic description we 
are using above would dominate over the non-equilibrium zero-sound modes. However, e.g. 
the modeling of ultra-high-density neutron stars is usually performed this hydrodynamic 
approximation . 

Having seen the level of complexity allowed even for the simplest systems one can 
ask what are the natural structures of the function G {X)? The simplest option for the 
diffusivity k is just to be constant. This corresponds to 



Another option would be to declare the diffusion coefficient D to be constant. This would 
correspond to G's and K's being connected as 



where we have omitted a constant of integration. In this case, G is just proportional to 
the charge density in equilibrium. However, as the simple physical example with fermions 
teaches us. Nature may not enjoy such structural simplicity at all. 

6 Discussion and Future Directions 

In this paper, we used a hydrodynamical language to describe the dynamics of theories with 
kinetic gravity braiding. This fluid picture turns out to be extremely useful: it provides 
an intuitively clear physical meaning to otherwise obscure combinations of derivatives of 
the field and the Lagrangian. This drastically simplifies notation allowing a much better 
understanding of the system. As is well known, finding the correct variables in many cases 
provides an easy path to the solution. Moreover, this notation clarifies the information 
encoded in the free functions K {(p, X) and G ((/>, X), making it possible to construct these 
theories based on physical principles and not solely from a naive naturalness. 

The key feature of our fiuid, is the dependence of the charge density and, corres- 
pondingly, the energy density on the expansion 6. This allows for the presence of bulk 
viscosity-like effects and diffusion without any dissipation. This novel property, to our best 
knowledge, had not been considered in the literature of hydrodynamics hitherto. 

We would like to mention that our identification of shift charges with particles and, 
therefore, of the proper-time gradient of the scalar field with the chemical potential is 
not actually unique. In fact, we could also have assumed that these charges correspond to 
entropy. From the first law of thermodynamics, it follows that this alternative identification 

^®We thank Gregory Gabadadze, Massimo Porrati and Alberto Nicolis for discussing this with us. 



G (X) = —\nX, for k. = const . 
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could be obtained by the exchange n — )■ s and m ^ T, where T is the temperature. We 
find that the identification presented in the paper is more satisfactory. In the case with 
entropy, even for a perfect fluid EMT, a non-shift symmetric Lagrangian would imply that 
entropy is not conserved. Moreover, even in the shift-symmetric case, our diffusion fiux 
would correspond to a heat flux, but one which would nonetheless conserve entropy. 

The equation of motion for the scalar is a highly nonlinear partial differential equation 
of the Ampere-Monge type. This class of equations belongs to the frontiers of current 
research in mathematics. In fact, contrary to the quasilinear case, it is not even clear 
under which conditions the Cauchy problem is well posed. It may happen that this fluid 
picture could aid in developing the understanding of this problem. In the forthcoming 
paper, Ref. [77] we make a first step in this direction by discussing the high-frequency 
stability and acoustic geometry in these theories. 

As we have alluded to already, our hydrodynamical picture owes a lot to Ref. [63], 
where perfect fluids with vorticity and thermal effects are described through a Lagrangian. 
It would be interesting to investigate whether a non-pathological and meaningful gener- 
alisation of perfect fluids with vorticity and thermal properties is possible along the lines 
discussed in our paper. We believe that for this purpose the Lagrangian V {(p, m, 6) from 
Eq. (B.IO) could be rather suitable. In Ref. [63] the author had already conjectured: "It 
may also be possible to extend this work to viscous fluids...". In fact, the current paper 
can be considered as a first step in this direction. Moreover, one could consider similar 
generalisations in the framework of Ref. [84] ; for a recent development in this direction see 
e.g. Refs [85, 86]. 

In a rather natural extension, one could employ the formalism developed in this paper 
to consider the hydrodynamical properties of the more general galileons and theories such 
as those of Ref. [53]. There is, however, a potential difficulty in interpretation which is 
related to a necessary appearance of curvature terms in the action, see Ref. [27, 87]. 

Further, one should not forget that boundary terms are in fact unavoidable for the 
proper formulation of the action principle in theories with higher derivatives, see e.g. [88, 
89]. It would be very interesting and important to find these terms and their possible fluid 
interpretation in theories with kinetic gravity braiding. 

Having in mind the higher-derivative structure of the theory, it is interesting to specu- 
late how one could change the fluid picture using constraints as in Ref. [90]. It seems that, 
in particular, one could obtain a kinetically braided dust — a braided fluid moving along 
geodesies. Another option would be to constrain the expansion in the action and obtain in 
that way a braided incompressible fluid with diffusion. 

Finally, the fluid presented here appears to naturally act so as to screen the grav- 
itational effects of matter when it is in an equilibrium conflguration. This is another 
manifestation of the monitoring effect that we discussed for the cosmological solutions of 
these models in Ref. [56]. Such phenomenology would clearly be important in cosmology 
during non-linear structure formation. Could the equilibrium conflgurations which we have 
described here provide a model for dark-matter haloes which significantly deviates from 
the standard CDM paradigm? 

To conclude, we are in these days witnessing a revival of interest and a rather exciting 
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series of developments in our understanding of non-canonical field theories and hydro- 
dynamics. The bestiary of models used in cosmological models has grown and we are in 
the process of trying to understand the implications. Eventually, these models will need to 
stand up to observational tests, whether in cosmology or possibly even in some condensed- 
matter systems, and we are in the process of trying to understand whether they could have 
any completely new signatures for which to hunt. 
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A Choice of Frames 

The choice of defined in Eq. (3.1) as a velocity for the frame is special in that it is 
explicitly vorticity free as a result of this velocity's being parallel to the gradient of cj). This 
means that we are permitted to perform a foliation of the space-time using (j) as an internal 
clock and the surfaces S<^: (j) = const as our spatial hyper surfaces. This particular foliation 
is then a natural candidate for being a Cauchy surface. In the follow-up work Ref. [77], 
we discuss the requirements on the possible general fluid configurations which ensure that 
this particular choice of spatial hypersurfaces indeed provides a surface on which initial 
values can be supplied in the usual unconstrained fashion and one can honestly describe 
the scalar field as a fluid. 

Using a general frame moving with some other velocity we could have defined the 
pressure, energy density and energy flow using equations (3.30), (3.32) and (3.33) with the 
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substitution — t- U^. In these variables the EMT would be 

= euU^U^- X-l^ Vu + q^^U^ + + vr^^ , (A.l) 

where 

TT - ( \ U \ U _\ \ U \ U \ rpaf} (A 2^ 

is the shear-stress tensor. Note that in principle we could perform this decomposition for 
an arbitrary EMT. Generically, none of the shear, rotation or vorticity will vanish in such 
a frame and the shear-stress tensor will contain terms with V and, correspondingly, rh. 
Hence, for general initial data, tt^i, will not vanish. In particular, this is the case for the 
two most popular frames used in the analysis of standard theories of imperfect fluids: the 
Eckart frame, where the frame moves with the particle flow (so that = / \/j\J\), and 
the Landau-Lifshitz frame, which moves together with the energy so that T^U^^^ = eLL^^LL- 
It could well be that these frames may not even exist for some otherwise reasonable initial 
data. 

One should note, however, that up to (and including) the first order in k, these two 
frames coincide, since the actual heat flux is absent, see [82, Eq. (2), p. 312]. Moreover, 
up to this order, the energy-momentum tensor in this Landau-Lifshitz frame takes the 
perfect-fluid form. And indeed, up to and including the first order in k we are then dealing 
with a perfect fluid. However, once higher-order corrections are included, we see that the 
price paid for this frame change is that the fluid flow is now no longer vorticity free, with 
corrections also coming at second order in k. We show this explicitly below. 

Let us for simplicity concentrate on the shift-symmetric case. In the Eckart frame 
(suppose it exists so that (3.19) holds), the shift-charge density is 

riE = = n + (^2) =p^^ + Ke + (^2) . (A.3) 

The 4-velocity of this frame is 

U^ = Jt'/nE = u'' + 0{k) . (A.4) 

Now let us calculate the kinematical rotation vector for the Eckart frame, see e.g. [91]: 

(f/^) ^ £"/37M (v,C/e/3) t/E7 = (J) /nl . (A.5) 

This kinematical rotation vector is related to the previously mentioned rotation tensor or 
twist in this frame w^j^ = 1/2 {^I^^V^Uev — -L[4tV"C/E/x) in the following way: 

{Ue) = [±'^^V^Ue(s) Ue^ = e^^^^^wI^Ue^ . (A.6) 
Further for uj^ (J) we have 

(J) = ^a/37A«v„ . (A.7) 

Some tedious but straightforward algebra yields the result that the vorticity in the Eckhart 
frame, in terms of the variables in our chosen local rest frame defined by the gradient of (/>, 

^""{Ue) = t^e'^f'^^So^Uf^a^ , (A.8) 

^E 
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with 

S^^^-V^6. (A.9) 

m 

Thus the kinematical rotation vector cj^ (Ue) is O (k^) and vanishes on equihbrium config- 
urations. In order to calculate it, one must solve the equation of motion (3.53). Nonethe- 
less, for general initial data, neither Vam nor Vq0 are parallel to a^, therefore u'^ (Ue) 
does not vanish. Nonvanishing vorticity implies (see Frobenius theorem) that the Eck- 
art frame is not suitable for the formulation of the Cauchy problem, because is not 
hypersurface-orthogonal. 

Using the variables we defined in section 3.7, we can also rewrite the full EMT 

T^u = (So + P)Ufj.Uy - Qf^i^P - Km{u^au + u^af,) + ©(k^) . (A.IO) 

We can now perform a diagonalization into the Landau-Lifshitz frame by boosting the 
velocity along a^, 

Ui^Lt, = Uf,-aaf, + 0{K^) (A.ll) 

with a = Km/{EQ + P). Rewriting Eq. (A.IO) using the new velocity, we obtain an EMT 
of perfect-fluid form, 

T^, = (So + P)ULL^^ULLu - g,.uP + 0{k^) . (A.12) 

However, at second order in k, deviations away from a perfect fluid reappear in this frame 
in the form of a non- vanishing anisotropic stress proportional to a^Ojy. 

The formulation of the theory in terms of the field (p provides a natural reference frame 
which is neither the Landau-Lifshitz frame nor the Eckart frame: the bulk velocity is 
neither the velocity of energy nor the velocity of particles. Contrary to the usual imperfect 
fluids, in this frame, the energy density and particle density contain first derivatives of this 
four-velocity. Moreover, the pressure contains the second time derivative of the scalar field. 
It is also in this reference frame where the theory has a cosmological solution and reduces 
to perfect-fiuid configurations. Finally, the shift-charge/particle-number current depends 
on initial data only and it seems that there are no general restrictions on the functions 
G {<j), X) and K (0, X) such that the current is timelike for all admissible initial data. 

B Action Without Second Time Derivatives 

Let us start from classical mechanics and consider the following one-dimensional version of 
the system represented by our action (2.1), 

S = J dt{K{q,q) + G{q,q)q) . (B.l) 

Now we can add a total derivative to this action without changing the equation of motion 
so that: 

qG {q, q) + ^/ {q, q) = q{G + f^) + 1,4 ■ (B.2) 
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We can eliminate the second derivative q if we take 

dvG{q,v) , (B.3) 

in which case the new, equivalent, Lagrange function without higher derivatives is 

Li {q, q)=K {q, q) - q J' dvG^, {q, v) . (B.4) 

Now let us see how one can eliminate the second time derivative for our kinetically 
braided scalar field (p. Let us use the Lagrangian (2.3) and write the action in the comoving 
frame in the form 

S^= I d^x / dr^^ (P - Km) , (B.5) 



where now g = g'^'^ _L with _L= det -L^u and 

g-- = g^'^d^rd.T . (B.6) 
Recalling the definition of the effective mass m Eq. (3.22) we have 



dr = — , so that g'^'^ = ^ = 1 . 



-,2 



Now we can add an arbitrary total derivative d yfV— -Lj /dr without changing physics so 
that 

5^1 = J^'^xJ drV^ (^P - Km + f^m + f— In + fmmj . (B.7) 

Therefore for 

rm 

f(cp,m)= dm'K{m',(j)) , (B.8) 



we eliminate m = (p from the action. Further, we notice that the expansion can be expressed 
as 

d = ^lnV^. (B.9) 
ar 

Thus we obtain an equivalent new Lagrangian V without second time derivatives in co- 
moving coordinates 

5^1 = J d^xV^ {P + mU + fO) = J d^x^V {(t>, m, 6) . (B.IO) 

Finally one can check that the difference between Lagrangians is indeed a total derivative: 

/m rm / rm \ 

dm'K4, + dm'K = V^A dm' k {m! ,4)) . (B.ll) 

Of course if we were to change the coordinates away from this natural LRF frame, V would 
again contain second time derivatives. 
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